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ABSTRACT 

We study how coupling to supergravity affects the phase structure of a sys- 
tem exhibiting dynamical supersymmetry breaking in a metastable vacuum. 
More precisely, we consider the Seiberg dual of SQCD coupled to supergravity 
at finite temperature. We show that the gravitational interactions decrease 
the critical temperature for the second order phase transition in the quark 
direction, that is also present in the global case. Furthermore, we find that, 
due to supergravity, a new second order phase transition occurs in the meson 
direction, whenever there is a nonvanishing constant term in the superpoten- 
tial. Notably, this phase transition is a necessary condition for the fields to 
roll, as the system cools down, towards the metastable susy breaking vacuum, 
because of the supergravity-induced shift of the metastable minimum away 
from zero meson vevs. Finally, we comment on the phase structure of the 
KKLT model with uplifting sector given by the Seiberg dual of SQCD. 
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1 Introduction 

Understanding how supersymmetry breaking occurs is a major problem on the road to 
connecting the underlying supersymmetric theory with the observed worlcfl The idea, 
that supersymmetry can be broken due to dynamical effects [2], has long been considered 
phenomenologically very promising, since it naturally leads to a large hierarchy between 
the Planck and the susy breaking scales. However, dynamical supersymmetry breaking 
has turned out to be quite difficult to implement in a supersymmetric gauge theory. The 

4 For a recent review of supersymmetry breaking see pQ . 
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reason is that only rather complicated examples [3] satisfy the strict conditions, necessary 
for the absence of a global supersymmetric vacuum. 

It was realized recently [3] that the situation changes dramatically, if one abandons 
the prejudice that the phenomenologically relevant vacuum has to be a global, and not 
just local, minimum of the effective potential. In this case, one can relax the requirement 
that the theory lacks a global supersymmetric vacuum and search for models with meta- 
stable, sufficiently long-lived, susy breaking vacua. From this new point of view, models 
with non-zero Witten index and without a conserved U(l) R-symmetry can be considered 
phenomenologically viable for supersymmetry breaking. The spectrum of susy breaking 
theories is then significantly enriched. In particular, as shown in [3], meta-stable dynam- 
ical supersymmetry breaking occurs even in Af = 1 SQCD with SU(N C ) gauge group 
and Nf massive fundamental flavors. This can be established by going to the Seibere, 
(magnetic) dual description of this theory, where supersymmetry is broken at tree-leveljj 
For convenience, we will call this the ISS model from now on. During the last year, many 
more examples of meta-stable dynamical susy breaking were found in various phenomeno- 
logically appealing settings [5]. Progress was also made on understanding the embedding 
of those field-theoretic models in string/M-theory [6]. 

However, once we consider phenomenology in a local, instead of a global, minimum of 
the zero-temperature effective potential, the following question arises. How natural is it for 
the high-temperature system, that is the early Universe, to end up in the metastable state 
after cooling down? To address this question, the recent works [7J[B] studied the ISS model 
at finite temperature. They found that the metastable vacuum is thermodynamically 
preferable compared to the supersymmetric global onesj^l Although their conclusions 
agree, their approaches are different and, in a sense, complimentary. In [TJ, they consider 
a path in field space, which extrapolates between the susy breaking vacuum and a global 
vacuum, and construct the effective potential along this path. Using this they show 
that, even if at high temperature the system starts at a susy vacuum, it will end up in 
the metastable one as it cools down. On the other hand, [8] studies in great detail the 
phase structure around the origin of field space, which is a local minimum of the nonzero 
temperature potential. They assume that at high temperature the quark and meson fields 
of the ISS model are localized near this point, which is reasonable since the number of light 
degrees of freedom at the origin is largest and hence this state maximizes the entropy. 
With this starting point, [S] investigates the phase structure of the free-energy as the 

5 We will review more details on that in the next Section. 
6 There are N c of them as we review in Section [5] 
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temperature decreases and finds out that there is a critical temperature, T®, for a second 
order phase transition in the quark direction (that is, the direction towards the metastable 
minimum). On the other hand, in the meson direction (i.e., the direction towards a global 
susy vacuum) they find that only a first order phase transition can occur, at temperature 
smaller than T®, and that it is quite suppressed by a high potential barrier J^| 

We follow the approach of |8J for the ISS model coupled to supergravity. This is 
certainly necessary for more realistic cosmological applications of the idea of dynamical 
supersymmetry breaking in a metastable vacuum. We will compute the one-loop effective 
potential at finite temperature by using the results of [2Tj, [22] for chiral multiplets coupled 
to supergravity. The nonrenormalizability of the latter theory is not an issue as it is 
supposed to be viewed as an effective low-energy description only, not as a fundamental 
theory. For more details on one-loop (albeit at T = 0) calculations in supergravity coupled 
with various matter multiplets, see [23J . However, we should note that the considerations 
of [2TJ [22] treat the zero-temperature classical supergravity contribution V to the T ^ 
effective potential as an effective potential itself. In other words, the T = loop 
corrections are viewed as already taken into account in the standard sugra potential 
Vq. So it might seem conceivable that these results may be affected (despite susy being 
broken) by regularization subtleties similar to those at T = 0. In the latter case, a 
regularization compatible with supersymmetry was developed in the last five references 
of [23] and was shown rather recently in [21] to sometimes have an impact on quantities 
of phenomenological interest, like the flavor-changing neutral currents. For the trivial 
Kahler potential, that we will need, this is not the case. Nevertheless, it is important, 
although going well beyond the scope of our paper, to address this issue in full generality 
at T ^ 0. Another remark is due. Every gravitational system exhibits instability under 
long- wave length gravitational perturbations [25J; this Jeans instability occurs also at 
finite temperature [26]. While it is certainly of great importance for structure formation 
in the early Universe, it is a subleading effect on cosmological scales on which the Universe 
is well-approximated by a homogeneous fluid. So we will limit ourselves to considering 
the leading effect, by using the formulae of [2"TII2"2"] for the effective potential, and will not 
address here the Jeans instability. 

We will show that the supergravity corrections decrease the critical temperature for 
a second order phase transition in the quark direction, T®, for any N c and Nf. While 
this is only a small quantitative difference with the rigid case, in the meson branch a 
significant qualitative difference can occur. The reason is that in the relevant field-space 

7 Reference [5] discusses in detail the suppression of this transition in a class of O'Raireataigh models. 
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region there are no contributions to the tree-level meson masses in the rigid limit and so 
the supergravity corrections are the leading ones. As a result, it turns out that, whenever 
the superpotential contains a nonvanishing constant piece Wo, there is a second order 
phase transition in the meson direction at a temperature Tjf, smaller than T®. However, 
this is not a phase transition towards any of the global supersymmetric minima, as it 
occurs at nonvanishing quark vevs. Notably though, our new phase transition is precisely 
what is needed for the system to roll towards the metastable vacuum since, due to the 
supergravity interactions, the latter is shifted away from the origin of the meson direction 
whenever Wo 7^ 0, as shown in [12] . 

Considering the ISS model plus supergravity is, in fact, the first step towards a full 
investigation of the phase structure of the KKLT scenario [10] with ISS uplifting sector 
at finite temperature. It was already argued in [ITj [T2] that metastable susy breaking 
provides a natural way of lifting the AdS KKLT vacuum to a de Sitter one, avoiding the 
problems encountered previously in the literature. Recall that the original proposal was 
to introduce anti-D3 branes, which break supersymmetry explicitly, whereas the later idea 
to use nonvanishing D-terms p2] turned out to be quite hard to realize [Hj. Studying the 
phase structure of the KKLT-ISS system is a big part of our motivation. However, in this 
case the computations become much more technically challenging. We make the initial 
step by showing that the origin of the ISS field space is no longer a local minimum of the 
temperature-dependent part of the effective potential. The shift of the high temperature 
minimum away from the ISS origin is related to the vev of the KKLT volume modulus. 

This paper is organized as follows. In Section [2] we review necessary background ma- 
terial about the ISS model. In Section [3] we compute the one-loop temperature-dependent 
contribution to the effective potential or, equivalently, the free energy of the ISS model 
coupled to supergravity. To achieve that, in Subsection 3.1 we derive the mass matrices 
for both quark and meson vevs nonzero, that are coming from the F-terms; in Subsection 
3.2 we take into account the D-terms. In Section 4 we expand the general results of Sec- 
tion 3 in terms of the small parameter Mp 1 , where Mp is the Planck mass. This allows 
us to read off the leading supergravity corrections to the rigid theory, considered in [5] . In 
Section 5 we compute the critical temperature for a second order phase transition in the 
quark direction to 0(Mp 2 ). In Section 6 we show that there is also a second order phase 
transition in the meson direction and estimate the critical temperature for it. In Section 
7 we consider the KKLT model with ISS uplifting sector and argue that the origin of the 
ISS field space is no longer a local minimum of the high temperature effective potential. 
The shift of the minimum away from this origin is determined by the vev of the KKLT 
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volume modulus. In Section 8 we discuss the implications of our results for the phase 
structure of the ISS model coupled to supergravity and for the end point of this system's 
evolution at low temperature. We also outline open problems. Finally, in Appendix A we 
give some useful formulae for mass matrices near the origin of field space and in Appendix 
B we show that no new supersymmetric minima appear in a small neighborhood of the 
origin in the ISS model coupled to supergravity, in the field directions of interest. 

2 ISS model 

It was argued in [15] that SQCD with SU(N C ) gauge group and SU(Nf) L x SU(Nf) R 
flavor symmetry has a dual (magnetic) description in terms of an SU(Nf — N c ) gauge 
theory coupled to certain matter fields. When the condition Nf < 3N c /2 is satisfied the 
magnetic theory is IR free. The matter content of the dual theory comprises two chiral 
superfields q and q, that are transforming in the (Nf , 1) and (l,Nf) representations of 
the flavor symmetry group respectively and are in the fundamental and antifundamental 
representations of SU(Nf — N c ) respectively (and so are called quarks), and a gauge- 
singlet chiral meson superfield $ in the (N f , Nf ) flavor group representation. Hence the 
index structure of the magnetic quark and meson fields is the following: g?, q J a and $| 
with i — 1, Nf and a = 1, ...,N , where N = Nf — N c . For easier comparison with the 
literature, we will also use the notation N m = N and N c = N e , implying in particular 
that Nf = N rn + N e . In terms of this notation the above condition for IR free dual theory 
is N e > 2N m . In the following we will only consider this case. 

The magnetic theory has the following tree-level superpotential: 

W = /zTrg$<z-/i// 2 Tr$. (2.1) 

The second term breaks the flavor group to its diagonal subgroup and corresponds to a 
quark mass term in the microscopic theory (i.e., the original SU(N C ) gauge theory). The 
Kahler potential is the canonical one: 

K = Trg f g + Tr<fg + Tr$ f $. (2.2) 

The magnetic description can be used to prove the existence of a metastable vacuum, 
which breaks supersymmetry at tree- level [3]. Indeed, it is immediate to see that the 
F-term condition 

= h (qfqi - ^51) = (2.3) 
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cannot be satisfied as the matrix qfqi has at most rank N m while 8\ has rank Nf. The 
moduli space of metastable vacua can be parameterized as: 




? T =I Z I > $ = I I " I > (2-4) 





where is an (iV/ — N) x (JVy — N) matrix while Q and Q are N x N matrices satisfying 
the condition QQ = yuI^xTV- The point of maximum global symmetry is at 

(qi) = (^)=^n, (g 2 ) = (g 2 ) = o, ($> = o, (2.5) 

where we have denoted: q T = ((71,(72) with qi and (72 being N x N and iV x (Ay — N) 
matrices respectively. It will also be useful for the future to introduce the following 
notation for the generic components of the Nf x Nf matrix <£>: 



$ = rU V12 , (2.6) 




where <pn is N m x N m , 12 is N m x N e , <f) 2 i is N e x N m and finally ^22 is an N e x N e 
matrix. 

The value of the scalar potential in each metastable minimum in (12.41) is: 

V min = (N f -N)h 2 fi 4 . (2.7) 

Usually, when supersymmetry is spontaneously broken, the moduli space of classical vacua 
is not protected against quantum corrections. As a result, the quantum moduli space is 
typically smaller and one may wonder whether any of the metastable vacua survive in it. 
In this regard, it was shown in [1] that the classically flat directions around the maximally 
symmetric vacuum (I2.5P acquire positive masses at one-loop through the supersymmetric 
Coleman- Weinberg potential |16j . This metastable minimum is therefore tachionic-free 
and from now on we will always mean (12. 5p . when we talk about a supersymmetry breaking 



vacuum. 

In addition to the perturbative corrections that we just discussed, there are also non- 
perturbative ones. Namely, gaugino condensation in the magnetic gauge group SU(N) 
induces the Affleck-Dine-Seiberg superpotential [17] : 

/ Hpt <f> \ w 
WAn S = N , (2.8) 



f\N f -3N^ 

where A is the UV cutoff of the magnetic theory, i.e. the scale above which the magnetic 
description becomes strongly coupled and hence not well-defined. Adding this dynamically 
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generated contribution to the classical ISS superpotential leads to N c super symmetric 
vacua, characterized by nonvanishing meson vevs£ 

3N-Nf 

h(<S>y ) = f ie—6j, (q) = (f) = 0, (2.9) 

in agreement with the Witten index [18J of the microscopic theory. The metastable 
vacuum can be made long-lived by taking e parametrically small as in that case the 
tunnelling to the supersymmetric vacuum is strongly suppressed]^] 

Since the ADS superpotential is suppressed by powers of the UV cutoff, for small 
meson fields it is completely negligible compared to the tree-level one, (12.11) . So in the fol- 
lowing we will drop Wads from our considerations, as we will study the finite temperature 
effective potential only in a neighborhood of the origin of field space. 



3 One-loop effective potential at nonzero T 

In the present section we compute the one-loop effective potential at finite temperature 
for the ISS model coupled to supergravity. Its analysis in subsequent sections will enable 
us to deduce the phase structure of this theory near the origin. 

Let us start by recalling some generalities about the path-integral derivation of the 
effective potential in a theory with a set of fields {x 1 }- An essential step in that is to 
shift x 1 by a constant background x 1 . Equivalently, we expand the Lagrangian around 
a nonzero background, {x 1 }, for the fields. Using this expansion, one can derive with 
functional methods a formula for the effective potential. The original derivation of [19] was 
only for zero-temperature renormalizable field theory. The same kind of considerations 
apply also for finite temperature and up to one- loop give [2D] : 

Veffix) = V tree (x) + V {1 \x) + V± l \x) , (3.1) 

where V tree is the classical potential, Vq is the zero temperature one-loop contribution, 
encoded in the Coleman- Weinberg formula, and finally the temperature-dependent cor- 
rection is: 

V?\x) = ~- (g B + ^f) + ^ [TrM s 2 (x) + 3TrM v 2 (x) + TWW?(*)] + 0(T) . (3.2) 
8 We denote by e the quantity fiA^ 1 . 

9 Recall that Nf > 3N and hence for very small e the meson vev, ($), in (|2.9p becomes very large. 
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For convenience, from now on we will denote this last expression simply by Vt- Here on 
and gp are the total numbers of bosonic and fermionic degrees of freedom respectively 10 !; 
TrM 2 , TrM 2 and TrM 2 are the coefficients of the quadratic terms of scalar, vector and 
fermiorj^] fields computed from the shifted classical potential or, in other words, the mass 
matrices of those fields in the classical background {x 1 }- The expansion (13. 2p is valid in 
the high temperature regime, more precisely when all masses are much smaller than the 
energy scale set by the temperature. The above result for the one-loop effective potential 
at finite temperature was shown in [21] to also hold for coupling to supergravity. 

We turn now to computing TrM 2 , TrM 2 and TrM 2 for our case. The classical back- 
ground x around which we will be expanding is: 

<9l> = (Qx) = Q ^N m xN m , (<Pll) = <pl ^N m xN m , (022) = <^2 ^N e xN e , (3.3) 

with zero vevs for all remaining fields and with Q, (fx, <p2 all being real. 



3.1 F-terms 

In this subsection we consider the contribution from the F-terms. The D-terms will be 
taken into account in the next one. For convenience, from now on we denote collectively 
all components of the fields q, q and $ simply by x 1 ■ 

3.1.1 Preliminaries 



Recall that the classical F-term supergravity potential isl 

V = e K {K n D 1 WD 3 W - 3|H/| 2 } , (3.4) 

where J, J run over all scalar fields in the theory, K IJ is the inverse of the Kahler metric 
and the Kahler covariant derivative is: 

DjW = diW + djKW . (3.5) 

The supergravity Lagrangian is invariant under Kahler transformations: 

Kix^X 1 ) - Kix^x^ + Fix^+Fix 1 ), 

Wix 1 ) - e-^Wix 1 ). (3.6) 



10 Note that this is different from the number of fields. For example, for Nb chiral superfields the 
number of scalar degrees of freedom is gs — 2Nb- In fact, below we will denote by Nb the number of 
complex scalars. 

11 In (|3.2[) the quantity TrM| is computed summing over Weyl fermions. 

12 In this section we set Mp = 1. The explicit dependance on the Planck mass will be reinserted later 
when needed. 
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One can exploit this invariance, by taking F{\ 1 ) — log V^(x 7 ) , i n order to show [27J that 
the scalar potential depends only on the combination 



G = K + \n\W\ 2 , (3.7) 
but not on W and K separately. In terms of this function, we can rewrite (13.41) as follows: 

V = e G (G J G/ - 3) (3.8) 

where 

ox 1 oxj 

This notation utilizes the fact that for us the Kahler potential is canonical i.e. Kjj = Sjj, 
see (E2D, and so xj = K jL x L '. 



For such a Kahler potential, the expressions for the scalar and fermionic mass matrices 
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are [2T] 

TrM s 2 = ( 2e G {(G IJ + G j G j )(Gjj + GjGj) + (N B - 1)G I G I - 2N B } ) (3.10) 

and 

TrMf = TiMl /2 + TrM 3 2 /2 = ( e G {{G IJ + G / G ,J )(G /J + GiGj) - 2} ) , (3.11) 

respectively. Here N B is the number of complex scalars. Using (I3.8p . it is easy to show 
that ( 13.101) follows from 

™- =2 ^pk- (3 ' 12) 

The derivation of the fermion mass-squared is a bit more involved since one has to disen- 
tangle the mixing between the gravitino and the goldstino in the supergravity Lagrangian. 
After having dealt with that, one finds TrM 2 ^ = —2e G . The anomalous negative sign 
in the gravitino contribution is due to the fact that in TrMjy 2 we have summed over the 
physical matter fermions and the goldstino. 

Before starting the actual computations, two remarks are in order. First, it may 
seem that it is more illuminating to perform the calculations as in [8], i.e. to compute 
separately the mass-squareds of every field and then add them. However, for us this 
becomes rather cumbersome, whereas the trace formulae above provide a very efficient way 
of handling things. And second, the formulation of the supergravity Lagrangian in terms 



13 We will consider a non-canonical Kahler potential in Section [71 when we address the KKLT set-up 
with ISS uplifting sector. 
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of the function G, (13.71) . appears to encounter a problem for vanishing superpotential, 
as W enters various terms in the denominator. That will be an issue for the D-terms 
in Subsection 13.21 and we will use there a more modern formulation that is equally valid 
for W 7^ and W = 0. Here we simply note that for the F-terms there is no problem, 
since the apparent negative powers of W, coming from derivatives of G, are cancelled by 
the positive powers from e G . (This will be made more explicit in Subsection 13.1.21 ) So 
the F-term results never contain division by zero. This is an important point as in later 
sections we will be interested in the effective potential at the origin of field space, where 
the ISS superpotential vanishes. 

3.1.2 Mass matrices 

We are finally ready to find the F-term mass matrices TrM 2 and TrM| for the ISS 
model coupled to supergravity. To do so more efficiently, we note that, instead of finding 
separate expressions for the various ingredients G ij Gjj 1 G I3 GjGj and G i G j GjGj which 
enter the mass formulae, it is computationally much more convenient to calculate the 
whole combination (G IJ + G I G J )(Gu + GiGj) for each particular choice of I and J. 
This avoids introducing a big number of terms that have to cancel at the end, as we now 
explain. From the definition of G, i.e. G = K + In W + In W, we have that 

Wj WrWj Wtt 

G, = K, + ^L, G, J = —L^ + 1 ±. (3.13) 

Therefore, in the expressions G^ G j j , G ^ G jG j and G^ G^ G jG j one seems to obtain 
many terms proportional to 1/W 2 W . Taking into account the \W\ 2 factor coming from 
e G , one is left with many terms ~ l/liyj 2 . However, it is clear that they have to cancel 
at the end, since the scalar potential (13.41) does not include any negative powers of \W\ 2 
and so if we were computing the masses of each field separately and adding them (as in 
[$]) we could not possibly obtain terms ~ 1 / 1 | 2 . This cancellation can be incorporated 
from the start by using (13.131) to write: 

Wtt KtWj + KtWj 

Gjj + GjGj = ~^f+ 1 w J 1 + KjKj (3.14) 

or equivalently 

Gjj + GjGj = ^L + G I K J + K I Gj- KjKj . (3.15) 

It is evident now that this expression does not contain any 1/W 2 terms. Incidentally, this 
also makes it obvious that, as expected, the expression e G {G IJ + G I G J )(Gu + GiGj) 
does not contain any powers of W in the denominator. 
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To illustrate how much the use of (I3.14p or (I3.15P simplifies the computation, let us 
look for instance at the following term: 



To make use of (I3.15p . we note that: 



(3.16) 



W A 



Therefore, one immediately finds: 



°) C^n) = <Pl %N m xN m , 

h 



u l 1 



( G (<t>22)f) 



'r2 



W 



$1 ■ 



(3.17) 



^011022^ 



HQ 2 - ^ 
Wo 



<P2+ [<P2 



w 



n 2 



(SiS l k )(5W), (3.18) 



where we denoted by Wo the value of the ISS potential in the background (13. 3p . The last 
factor gives (<5?4)(^f) = = N mN e - So we obtain 

2 



¥1^2 + 



h(Q 2 - /i 2 ) hn 2 



Wo 



-^2 



-<Pi 



N m N e . 



(3.19) 



Similarly, it is very easy to compute: 



2 

m 1 



+ 2 



hQ 2 
Wo 



, h(Q 2 -^ 2 



Wo 
h 



+ 1)<P1 + 



h(Q 2 - 12 
Wo 



N„ 



¥2 + —{Vl¥2-^ 



N m N e , 



— I N m N e , 
Wo ' 



2h 



(R^) = (R qiqi ) = Q 4 1 + — Vl N\ 



Wo 



I 2 

m ) 



(R 919-1 ) 
(R q2q2 ) 



Wo 

Wo 



A^ + Q 4 1 + — fi + 



2/i(^ig 2 / 2/iy? 



Wo 



Wo 



N N 



N m 
(3.20) 
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For all remaining pairs (R x x ) = 0. 

The last ingredient in (13.101) . that we need to compute, is (G I Gi). The only nonvan- 
ishing components are for / = </>n, 2 2, qi, q~i and we find: 

{G , Gj) = ^ + ^1 ) \ m + (V 2 . 5g) 2 + 2Q2 + t ) " Wm . ,3. 21 ) 

One can notice that all expressions above depend only on Q 2 , not on Q alone. Hence 
the temperature dependent part of the one-loop effective potential as a function of Q is 
of the fornix 

V T = e Q2 (AQ 6 + BQ A + CQ 2 + D) (3.22) 

for any values of the meson vevs ip% and ip2- Therefore Q = is always an extremum. In 
the meson directions things are not so apparent, as there are odd powers of (pi and <p2- 
However, one can see that all of them multiply either a power of Q 2 or the first power of 
Wo- Since Wo is linear in the meson vevs, for Q = the dependence of Vr on ifi and <p2 
is, in fact, at least quadratic (by that we also mean mixed terms, i.e. with ^1^2) • We 
will see in Section fl~Tl that the point (Q, ipi,(p 2 ) = (0,0,0) is a local minimum of Vr, as 
was also the case for vanishing supergravity interactions Before that, however, let 

us first consider the D-term contributions to the mass matrices of the various fields. 



3.2 D-term masses 

The D-terms for Super Yang-Mills coupled to supergravity were derived for the first time 
in [27]. However that formulation, entirely in terms of the function G of eq. (13. 7p . is 
not convenient for our purposes since the D-terms have a singular dependance on the 
superpotential. For example, denoting by g the gauge coupling constant and by T a the 
generators of the gauge group, the D-term scalar potential was found to be 

V D = ^D a D a , (3.23) 

with 

D a = 9 Gj T a 5 X J = \g T a k J ^, (3.24) 

where in the last equality we have substituted G 1 = K 1 + W 1 /W. Clearly, V D is not 
well-defined when the superpotential W vanishes. The same problem, i.e. division by W, 

14 We drop from now on the piece that is ~ T 4 , as it does not depend on the vevs of the fields and so 
it does not contribute to the derivatives of Vr, which are the quantities that will be of interest for us. 
15 This is no longer true if one includes the KKLT sector, as we argue in Section [7] 
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appears also in the fermionic mass terms. So if we adopt the formulation given in 
TrMj seems to diverge for W = 0. That is problematic since we would like to study the 
effective potential at the origin of field space, where the ISS superpotential vanishes. 

It has been noted long ago that the above formulation is not suitable in the presence 
of a vanishing superpotential and the latter case has to be studied separately, without 
the use of the function G = K + In \ W\ 2 . However, a careful study of the supergravity 
Lagrangian that is valid both for W 7^ and W = was performed only recently (to the 
best of our knowledge) in [28]. What is relevant for us is that Vd can be written as 

V D = ^ I a d I K-^r a )\ (3.25) 

where are the gauge transformations of the scalar fields, i.e. SaX 1 = £a({x J })> an d r a 
are functions determined by the gauge variations of the superpotential: 

S a W = CdiW = -3r a W , (3.26) 

where the second equality is required for gauge invariance of the action. These functions 
also characterize the gauge-non-invariance of the Kahler potential: 

6 a K(x,x) = Hr a ( X )+f a (x)). (3.27) 

In the case of non-vanishing superpotential one can use (I3.26P to express r a in terms 
of W and d{W . Substituting the result in (13.251) . one finds (I3.23p - fl3.24p upon using 
£^ = iT a jX J - This, indeed, shows that the formulation of [2E] reduces to the one in [27] 
when W + 0. 

Since the ISS superpotential is gauge invariant, we have r a = and 

2 ^m-l 2 

Vd = Y S (M4T a qi-qiT a q\)) (3.28) 
0=1 

as in [8]. We can then borrow the results, found in the global supersymmetry case, to 
obtain a 4g 2 Q 2 (N% l — 1) contribution to TrM 2 . Also, the vector boson mass is the same 
as in [8] and so it gives a 4g 2 Q 2 (N^ — 1) contribution to TrM^. 

Let us now come to the fermionic sector. The mass matrices are |28|: 



M u = v i V J M (329) 
M Ia = -i 

M a/3 = -]djf a/3 K IJ Mj, (3.31) 



d!V a - ^(Refy^VpdJ^] (3.30) 
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where f a p are the gauge kinetic functions, the action of the covariant derivative V 1 on 
M = e K ' 2 W is 

V 1 = d T + -d I K, (3.32) 

M 1 = V T M and finally 

V a = i&d I K. (3.33) 

These expressions are clearly well defined and non singular even when W = 0, unlike the 
analogous formulae in [27] . In our case, obviously the index I runs now only over the 
quark fields. For a flat Kahler potential, one can easily verify from (I3.29p that: 

M IJ M U = e G (G IJ + G / G' J ) (G u + Gfij) , (3.34) 

so that we recover correctly the contribution from the matter fermions and the goldstino 
in (13.111) . Since for us f a p = 1/g 2 = const for all a and /3, M a p = while the contribution 
to TrMj from the mixing of gaugino and hyperini can be found from (13.301) using V a = 
q[T a qi — qiT a q{. It reads 

(2M Ia Mi) = Sg 2 Q 2 (Nl - 1) . (3.35) 

In the supergravity Lagrangian there is one more mixing between fermions, which could 
potentially add a term to TrM|, namely the mixing between the gravitino and gaugino. 
In [28j it is of the form %l)V a \ a . However, in our case (V a ) = and so this mixing does 
not contribute. To recapitulate, the D-terms give exactly the same contribution as in the 
rigid case considered in [8]. 



4 Expansion in M p l and rigid susy limit 

In Section [3l we computed all ingredients of the finite temperature one-loop effective 
potential for the ISS model coupled to supergravity. Now we will make connection with 
the globally supersymmetric theory by inserting in the relevant formulae the explicit 
dependance on the Planck mass, Mp, and expanding in powers of Mp 1 . For later use, we 
will extract the leading supergravity corrections to the rigid results. 
Let us start with the tree-level supergravity potential 

V = Mp e G (M P G I G I - 3) + l -D a D a (4.1) 

with 

K \W\ 2 

G = ——77 + log . (4.2) 

Mp S Mp K J 
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As discussed in the previous section, the D-term contribution is the same as in the rigid 
limit. It is easy to see that the expansion of (14. II) gives: 

V = W T Wi + Vd" + —2 {KW I W I + 2Re(K I W I W) - 3\W\ 2 ) + 0(Mp A ) , (4.3) 

where obviously the first two terms are the standard global susy result. For future use, we 
note that taking / = n , 2 2 in the last equation gives, to leading order in the supergravity 
corrections, the following contribution to the classical F-term potential for the quarks: 



V F D (h 2 Tr q 1 q\q\q 1 - h 2 // 2 (Trqm + Trgjgj) + h 2 ^(N rn + jV e )) (l + 



Trgig} + Trgiqj 
Ml 



(4-4) 

Inserting the Mp dependence in the thermal one-loop potential Vr yields 16 !: 

V T = ^M 2 (e G (^M P J2 R U + 2(iV B - ^M^G'Gj - 2{2N B + 1) j ) + (V D ) , (4.5) 

where we have denoted by R IJ the quantity {G IJ + G I G J )(Gu + GjGj), as before. To 
obtain the explicit powers of M P in R IJ , we note that due to (14. 2p equation (I3.15P becomes: 



W u , GjKj , GjKj KjK, 



G IJ + G lGj = ^ + ^ + ^-^. (4.6) 
Now, expanding (14.51) we find: 

V T = ^((3W IJ W U + ^W IJ W LJ + ^Re[^(W I K J + K*W J ) 

+ 2 {NB M j 1) W I W I )) + (V D } + 0(M p *). (4.7) 

Together, equations ( 14. 31) and (14 .7\\ give the general expression for the zeroth and first 
orders in the Mp 1 expansion of the one-loop effective potential at finite temperature. 

Let us now apply the above formulae for the background (13.31) . but with vanishing 
meson vevs. To see what (14 .71) leads to, note that the only contributions to the first 
two terms come from (\W qi(f>11 \ 2 ) = (\Wq l4>11 \ 2 ) = h 2 Q 2 N^ and (\W q2(hl \ 2 ) = {\W q2 ^ 2 \ 2 ) = 
h 2 Q 2 N m N e and that the only nonzero (W 1 ) are (W* 11 ) = h(Q 2 - /i 2 ) and (W^ 22 ) = -hfi 2 . 
Hence we obtain: 



Vt = ^h 2 Q 2 (N 2 m + N m N e )(l + 2^)+T 2 Q 2 g 2 (Ni-l) (4.8) 
T 2 / I 



+ [h 2 Q 2 (Q 2 - fi 2 )N m + -{N B - l)[h 2 (Q 2 - fi 2 ) 2 N m + h 2 ^N e ]j + 0{Mp A ). 

16 Recall that, as we mentioned in Footnote [TU we drop for brevity the constant ~ T A contribution to 
effective potential in (|3.2[) , as it does not affect our considerations. 
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Taking Mp — > oo, we find the global supersymmetry result for wf^l 
4.1 Local minimum at the origin 

At high temperature, the temperature-dependent contribution Vp completely dominates 
the effective potential and so the minima of V e ff are given by the minima of Vt- Let us 
now address the question whether the origin of field space is a local minimum of Vt- 

In principle, this could be a complicated problem, as we have to consider a function 
of three variables, Vr(Q, fi, y^)- However, things are enormously simplified by the fact 
that, as we noted at the end of Subsection 13.1.21 Q = is an extremum for any (p% 
and if2- Since it is also a local minimum in the rigid limit, it will clearly remain such 
when taking into account the subleading supergravity corrections in our case. So we are 
left with investigating a function of two variables, tp\ and (f2- The presence of terms 
linear in any of them could, potentially, shift the position of the minimum away from the 
point (<fi,(p 2 ) = (0,0)o However, such terms in the ISS model coupled to supergravity 
appear only multiplied by powers of Q 2 , as we noted below eq. (I3.22p . So, given that 
(ifi, if2) = (0, 0) is a local minimum in the global supersymmetry case, it will remain such 
after coupling to supergravity, to all orders in the 1/Mp expansion. 

It is, nevertheless, instructive to write down explicitly the expression for Vt to leading 
order in the 1/Mp corrections: 

V T (Q,<pu<P2) = ^h 2 Q\Nl + N e N m )+T 2 Q 2 g\Nl-l) + ^h\^\Nl + ^ 2 N m N e ) 



+ 



T 2 |~.o„o,_ o „o o, (Nb — 1) 



M 2 



h 2 Q 2 {2 v \ + Q 2 - ^)N m + v B u > {h 2 {{Q 2 - /i 2 ) 2 + 2tf)N n 



! h 2 fi 4 N e } + (2Q 2 N m + V \N m + V \N e ) ( h 2 Q 2 (N 2 m + N m N e 



2 



+ h 2 ( V lN 2 m + V 2 2 N m N e 



1 



M P 



+ ^(774)- (4-9) 



Clearly, this is consistent with (14.81) . We have collected the terms that survive in the 
Mp — > 00 limit on the first line. Obviously, the origin of field space (Q, (fx, ^2) = (0, 0, 0) 
is a minimum. One can also notice that to this order Vt is a function of the form 
Vt(Q 2 , (pi, y? 2 ), i.e. it does not depend on odd powers of any of the vevs. One can verify 



17 Our result is in agreement with that of [5], upon correcting a typo there. Namely, they have over- 
counted by a factor of two the number of Weyl fermions (in their notation) ^ff x and (^f + \&f ) / \/2. 
18 This is what happens for the KKLT-ISS model, as we will see in Section [7] 
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that the terms linear in (p and multiplying powers of Q 2 , that we mentioned above, start 
appearing at 0(1/ Mp), while terms linear in (px(p% appear at 0(1/Mp). 

It will be useful for the next sections to extract from (14. 91) a couple of special cases. 
One case is the second derivative in the quark direction only: 



d 2 V T 



=0,ipi=0,ip2=0 



T>[h\Nl+N m N e )+2g\Nl-l)Y 



T 2 1 



M 2 3 



-h 2 ^N m (N B +5). (4.10) 



The other interesting case is the potential in the meson directions only: 



V T (0, ipi, if 2) 



+ 



H 

4 

1 1 



h 2 



h 2 (rtN 2 m + ip 2 N m N e ) + w (ip\N m + <f%N e yN t 



M 2 3 



-h 2 ^(N B -l)(N m + N e 



(4.11) 



Finally, let us note that if we want to study other minima of the potential, which are 
away from the origin of field space (Q, (fi, cp 2 ) = (0, 0, 0), we would have to include in our 
considerations the non-perturbative superpotential Wads, see ec i- 



5 Phase transition in quark direction 

As we saw in the previous section, at high temperature the origin of field space is a local 
minimum of the effective potential. Lowering the temperature, one reaches a point at 
which this minimum becomes unstable and the fields start evolving towards new vacua. 
As recalled in Section [21 at zero temperature the ISS model possess supersymmetric 
vacua at non-zero meson vev and a metastable vacuum in the quark branch. Adding 
the supergravity interactions results in a slight shift of the positions of those minima. In 
particular, the metastable vacuum shifts to small vevs for some of the meson directions 
[T2j . Hence, in order to end up in it, the system must undergo a second order phase 
transition in those same meson directions, as well as in the quark ones. In the next section 
we will see that this is indeed the case. In the present section, we will find the critical 
temperature, T®, for the onset of a second order phase transition towards nonvanishing 
quark vevs, while keeping all meson vevs at zero. 

Before turning to the supergravity corrections to T®, let us recall how to compute 
the critical temperature in a generic field theory [2U]. Suppose that we have a theory 
with a set of fields {x 1 }- To find the effective potential, one has to shift x 1 by constant 
background fields x} ' , as we reviewed in Section [3) The effective potential is a function of 
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X 1 only and we consider the case when it is of the form V e //(x 2 )o The location of the 
minima is then determined by: 

Clearly, x 1 = satisfies this condition. Other minima, at % / 7^ 0, can only occur when 
dVeff/dx 2 = 0- The critical temperature, at which rolling towards such minima begins, 
can be found by requiring [20J: 



dV„ 



eff 



dV 



dV T 

+ 



= 0, (5.2) 



dx 2 x=o dx 2 x=o dx 2 

where we have split the one-loop effective potential into zero-temperature and temperature- 
dependant contributions, Vo and Vr respectively. This is equivalent tec! 

dx 2 lx=o 2 

where m 2 is the unshifted tree-leve|^ mass-squared of the field, whose nonzero vev char- 
acterizes the new vacuum. This last equation can only be solved if m 2 < 0. In the global 
supersymmetry case, among the quark fields, which have nonzero vevs at the metastable 
minimum, the only ones with negative tree-level mass-squared are the components of 
Re( qi +ql)/V2 = Req+ [8]. Since for small vevs the gravitational corrections are sub- 
leading, we are guaranteed that this will be the case for us as well. 



19 This will be the case for us. However, even in principle this assumption is not a restriction but merely 
a simplification, as we explain in a subsequent footnote. 

20 The assumption that V e ff depends only on x 2 , but not on \ alone, guarantees that \ — is an 
extremum of the effective potential and thus simplifies the computations technically. In the general case, 
instead of (|5.3|) one has to solve, symbolically, the following system: V^j(\ Cl T c ) = , V"^(x c , T c ) = in 
order to find the critical temperature T c , at which the phase transition occurs, and the vev \ c , at which 
the relevant extremum of V e ff is situated. Clearly, here we have denoted by ' and " first and second 
derivatives w.r.t. to x, respectively. 

21 In principle, m 2 gets also contributions from the Coleman- Weinberg potential. However, in our case 
the origin of field space is not a local minimum of the zero temperature potential and so perturbation 
theory around it does not make sense. Hence, for us, m 2 is purely classical. In fact, to be more precise, 
we should note that at the origin the Coleman- Weinberg potential becomes imaginary. As shown in 
[29] , the imaginary part of the effective potential encodes the decay rate of a system with a perturbative 
instability. This decay process leads to the breaking up of the initial homogeneous field configuration into 
various domains, in each of which the fields are rolling toward a different classical solution, and should 
not be confused with a non-perturbatively induced tunneling between different minima of the potential. 
We will not dwell on that further in the present paper. 



18 



5.1 Mass matrix diagonalization 

We will now compute the tree-level squared masses of the fields Req + , that are necessary 
for finding the critical temperature in the quark-direction. Before considering the super- 
gravity corrections, it will be useful to give the derivation of these masses in the global 
supersymmetry limit. For more generality, we keep Q ^ in this latter case, although we 
will take Q = when we turn to the supergravity contributions to m 2 in (15. 31) . 

From the tree-level scalar potential of the global theory, V = K IJ OiWdjW , one finds 

[8]: 

h 2 Q 2 (\qi\ 2 + \qi\ 2 ) + h 2 Q 2 Tr( qi q\ + gfe) + h 2 (Q 2 - fi 2 )(Trq iqi + h.c.) . (5.4) 

When we turn on gravitational interactions we will have a similar expression but, generi- 
cally, with different coefficients. So it is of benefit to consider the general expression: 

A(\<li\ 2 + \Qi\ 2 ) + B(Ti qi q\ + Trgjgx) + C(Trg l9l + h.c.) (5.5) 
with arbitrary A, B and C . It can be diagonalized easily by introducing the combinations 

9i + ?* Qi - q* /, a s 

^ = -^-. = — ■ < 5 - 6 > 

Substituting the inverse transformation, 

q+ + q- , ~ t q+ ~ 9- ^ 

q 1 = and q 1 = , (5.7) 

in (15.51) . we find: 

A{\q + \ 2 + |g_| 2 ) + B{\q + \ 2 - |g_| 2 ) + C \ Tr(g 2 - q 2 _ + $ - ql 2 ) . (5.8) 

Finally, we decompose q± = Re(q±) + ilm{q±) and obtain the following mass terms: 

(A+5+C')(i?e(g + )) 2 +(A+fi-C')(/m(g + )) 2 +(A-5-C')(i?e(g_)) 2 +(A-fi+C , )(/m(g_)) 2 . 

(5.9) 

Reading off the values of A, B and C from ( 15. 4ft . we obtain from ( 15. 9ft : 



,2 !,2/o^2 ..2 



Req+ 



h 2 (3Q 2 -^ 2 ) , m 2 mg+ =/i 2 (Q 2 + /U 2 



L,_ = W-Q 2 ) , ^L 9 _ = ^(Q 2 -^ 2 )- (s.io) 



We see that, as already mentioned above, only the components Req + of the field q + have 
negative mass-squareds for zero shift Qc3 Therefore, only their masses should appear on 



22 Clearly, the negative mass-squared of Iraq- is irrelevant as in the metastable supersymmetry breaking 
vacuum = 0. 
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the right hand side of eq. (15.31) . We note that the masses for the fields Imq + and Imq- 
differ from those given in table 5 of [8]. Fortunately their typos cancel out in the total 
TrM s 2 . 

Let us now apply the result ( 15. 9ft for the diagonalization of the expression ( 15.51) to the 
case of interest for us. Namely, we consider the ISS model coupled to supergravity and 
we want to compute the following derivatives: 

d 2 t y, d 2 ^V, dlVetc, (5.11) 
nil gill qiqi 

which give the coefficients A, B and C. As before, V is the tree-level supergravity scalar 
potential. The only nonvanishing derivatives, upon setting Q = 0, arecf 

m 2 t = ml _ t = —prh 2 ii\N m + N e ) , ml, = m 2 t _t = -h 2 fi 2 . (5.12) 
qM nil Mp qiqi 1W1 

Hence, applying (15.91) . we find that the tree- level mass of each of the real fields in 
Req + is: 

<e q+ = + j^h 2 ^(N m + AQ . (5.13) 

The first term corresponds to the global supersymmetry result, as can be seen from ( 15.101) . 
while the second is a correction due to the supergravity interactions. 

Note that we did not need to include the D-terms in the computation of the unshifted 



masses. The reason is that, as we have seen in Section [3T2l the contribution of these terms 
to the mass-squared of the fields is always proportional to Q and so vanishes for zero shift. 

5.2 Critical temperature 

To find the critical temperature in the quark direction, recall that we shift the fields qif 
and q~i b by constant matrices QSf and Q5 3 b respectively, see eq. ( 13 .3p . In the basis of q + and 
g_, see equation (15.61) . the only fields which get shifted are the N m diagonal components 
of q + . From those, only the N m fields Reg + have negative tree-level mass-squared as we 
saw above. Therefore we can find the critical temperature fronj^ 

d 2 V r o 

Q=Q = -N m m 2 Req+ . (5.14) 



23 For more details see Appendix A. Note that (|5 . 1 2[) are exact expressions, i.e. to all orders in 1/Mp. 
However, since they happen to be at most of O(-^p-), they can also be read off from the part of the scalar 
potential (|4.4|) . which contains only the relevant leading order supergravity corrections. 

24 Clearly, this is equivalent to IEM as d 2 V T /dQ 2 = 2dV T /d(Q 2 ). 
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Recall that, to leading order in the supergravity corrections, we have (see (14.101) ): 

f^L = t2 [ a2 « + ^ + ^« - 1)] - jjg 7a * v y' +ii) ■ (5-i5) 

Together with 05.13(1 . this then leads to: 

h 2 fx 2 N m 1 h 2 p, 4 (N m + N e )N m 



( T cf 



h\Nl + N m N e ) + 2g\Nl - 1) M 2 /i 2 ^ + N m N e ) + 2^ 2 (iV^ - 1) 
J h^Nl(N B + 5) / J_\ 

Mp 3 [/i 2 (A^ + N m N e ) + 2g 2 (N^ — l)] 2 l^MpJ ' l ' J 



To see whether the order 1/Mj, correction increases or decreases the critical temper- 
ature, recall that iV/ > 3iV i.e. N e > 2N m . Let us write 

N e = 2N m +p, peZ + (5.17) 

and substitute this in the total numerator of the 0(jp-) terms in (15.161) : 

Numerator = -h 4 fi 4 Un^ + ^-N mP + 2p ~ 5 > | ^ - 2h 2 fx 4 g 2 (3N 2 n + pN m )(N 2 m - 1) . 



3 3 

(5.18) 

Clearly, for any N m and for any p > 1 every term in the above expression is negative 
definite. For p = 1 there is a positive contribution from the last term in the first bracket: 
h 4 fi 4 N^. However, it is outweighed by the first two terms in that bracket for any value 
of N m . So the conclusion is that the supergravity interactions cause to decrease 
compared to the rigid case for any N m and iV e . 



6 Phase transition in meson direction 

In the global supersymmetry second order phase transition is only possible in a 

field space direction with nonzero squark vevs [8]o However, in our case things may 
be quite different due to the supergravity-induced contributions to the tree-level meson 
masses. Note that, in the approximation of neglecting Wd yn , the meson fields have zero 
classical mass in the global limit (8] and so the supergravity corrections are the leading 
oneso Thus the possibility occurs that in the region of small vevs, for which neglecting 
the dynamically generated superpotential is well-justified, some of the meson directions 



25 Although a first order phase transition can still occur in the meson direction [5]. 
26 Recall that we do not have to take into account contributions from zero-temperature one-loop effects 
in the rigid theory as the origin of field space is not a local minimum of the zero-temperature potential. 
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may develop negative tree-level mass-squareds due to supergravity. We will see in the 
following when that happens. 

Let us address this issue in a slightly more general set-up as in [T2] , namely by adding 
to the ISS superpotential a constant piece. I.e., we consider W = Wo+Wjss, where Wjss 
is as in (12.11) and Wq = const. This is also a useful preparation for the KKLT-ISS model 
that we will discuss more in the next subsection. It will turn out that we need to compute 
the tree-level supergravity-induced meson masses not only at the origin of field space but 
also along the quark direction, i.e. for Q ^ and (pi,(f2 = 0. Ezl Also, we will have to keep 
track of terms of up to 0(1/Mp) to see the novel effect. 

Turning to the computation, from flA.61) we find: 



m (^n)i(0u)i 



m 2 



= 2h 2 Q 2 8j8f + ^ {5$ h 2 [(Q 2 - ^) 2 N m + ^N e ] 

+ 8j5f h 2 [(Q 2 - /i 2 )(3Q 2 + /i 2 ) + 4Q 4 iV m ] } 

+ ^ [~ 2W o + 2h 2 Q 2 N m ((Q 2 - ^) 2 N m + fN e )] 

+ 5j5f Ah 2 Q 4 N m [2(Q 2 - fi 2 ) + Q 2 N m ] 
1 



{515; h 2 [(Q 2 - ii 2 ) 2 N m + /i 4 iV e ] - h 2 ^5 n m 5l} 



p 

+ ^ &J} {-2W 2 + 2h 2 Q 2 N m [(Q 2 - ^ 2 ) 2 N m + /i 4 iV e ] } 



2 ; 2 i Q~ 4 /'"' . 4Q 4 iV m | 



W>ll)^2 2 )? ^ Ml Mi J * p ' 

Q 2 +^ 2 , 4(J 4 JV„ 



where clearly i,j,k,l = l,...,N m and m,n,p,q = l,...,N e . Note that VV appears only 
at order 1/Mp. All terms of the form ™\ xx ^ v m 22 </> 2 2 an d m 0n<^ 22 van i sri ( see Apendix 
A). To understand which mass matrix one has to diagonalize, let us recall that the meson 
fields are shifted as follows: (0n)| by ipi5j and (022)™ by <f2^, see eq. (13 .3p . Hence there 
are N m fields shifted by if\ (let us denote them by 0i) and N e fields shifted by tp 2 (let 
us denote them by ($2)- in the quark direction we could factor out an overall N m (since 
both gi and q~i have the same number of components) and diagonalize the remaining mass 
matrixo Now things are somewhat more complicated as there are different numbers of 



27 For reasons that will become clear below, we will be interested only in the range Q £ [0,/j,]. So we 

are still allowed to neglect Wd yn - 

28 This is why we ended up with m 2 = N m m 2 Req+ in (|5.14[) . where m 2 Req+ was the result of the diago- 
nalization. 
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and 02 fields. More precisely, from (16.10 we see that we have to diagonalize the expression: 



Ml = {2h 2 Q 2 N 2 m + ^(N m C l + N 2 m C 2 ) + ^(N m C 3 + N 2 m C,))<j )l ^ 



1 



( h 2 (Q 2 - ii 2 ) 2 N m N e + j^N e C 3 ) 0202 



where 



h 2 ii 2 N m N e 



Cl 

c 4 



Q 2 + fi 2 , 4Q 4 iV, 



M 2 



+ 



Mf, 



(0102 + 2 0l) 



(6.2) 



h 



(g 2 -/i 2 ) N m + fi 4 N e 
h 2 [(g 2 -// 2 )(3g 2 + // 2 ) + 4g 4 iv m ] 
-2iy 2 + 2g 2 iv: m c 1 , 

4/i 2 Q 4 iV m (2(Q 2 - / i 2 ) + g 2 iV m ) . 



(6.3) 



Notice that Wq enters the above formulae only at order 1/Mp. 

Before proceeding further, let us make an important remark about the position of the 
local minimum that is our starting point at high temperature. From (13.201) and (13.211) . one 
can see that for Wq ^ terms linear in (pi and (p% start appearing at order 1/Mp. This 
implies that the position of the minimum is shifted to some point (Q, if 1,^2) = (0, y*, ip£) 
with ipZ, 7^ 0. However, one can easily verify that both ip\ and are °f 0(1/Mp)\^_ 
Therefore, in all terms of (1A.1I) with an explicit 1/Mp dependence one can take V 9 !; = 
since we are working to order 1/Mp. Hence the only place the nonzero <p*, ip\ can make a 
difference in is the zeroth order terms. This translates to the zeroth order term in (1A.6j) . 
However, that term is independent of the meson vevs as the superpotential is linear in 
those and the indices /, J are along meson directions. Thus the results (I6.ip - (l6.3p are 
unchanged by the nonzero </?*, ip\. For convenience, in the following we will keep refering 
to the local minimum at (0, <p\, ip^) as 'the minimum at the origin of field space'. 

Let us now go back to (I6.2p . Since the coefficients of 0i0i and 0202 are not the same, 
we cannot diagonalize this expression immediately by using the results of Subsection 
5.11 However, it is still useful to change basis to the real components of the fields: 
0i = i?e0i + i/m0i and 02 = -Re02 + 2-Zm0 2 . Then M? acquires the form: 



Mj = J] (Mn x 2 + 2M12 + M 22 yt 



(6.4) 



J More precisely, one finds ipl 



_ in 2 W 1 
— hN m Ml 



O(j^) and & 



hN m Ml 



looked at the subleading orders, so we abstain from claiming that ip\ — ipi,. 



G(-L). We have not 
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where (xi,yi) = (Refa, Refa) and (£2,2/2) = (Imfo, Imfo). Thus the problem is reduced 
to the diagonalization of the 2x2 matrix M^ with a,b — 1,2 and so the eigenvalues are 
given by: 

2 TrM ± A/(TrM) 2 — 4detM 



(6.5) 



Before proceeding further with the Q 7^ considerations, let us first take a look at what 
happens at the point (0, cp*, (p^), which for us, as explained above, is the same as looking 
at the origin of field space. 

For (5 = 0, the expressions ( 16. II) and ( 16. 21) simplify significantly. Namely, we have: 

, h 2 ^ 



m, 



m 



L)i(0n)i lOr 



(022)^(<A22)?IOr 
m (0nH(022)?lor 



Mj, 

M 2 



5^f(N m + N e )-SlS, 



3xk 



2W? 



6 q 5? 

w m w p 1 



h 2 u 4 - 
Ml 1 p ' 



m 



(^22)^(^11)1 1 Or 



h 2 n 4 

M 2 m l ' 



(6.6) 



where 



Or 



M 



denotes evaluation at the origin, and therefore: 



I Or 



-2 iv m iv e ^4 ^V m 



2, ,4 



+ 



M 2 



2iy n 2 



h 2 n 



N m N, 



+ 



)• 



Applying (I6.5p . we find that the two eigenvalues are: 
2/iV 



M 



N N 

2 ly m 1 " e 



W l {N m + N e 



Mp 



m 



W i-(N m + N e 



M P 



(6.7) 



(6i 



So, clearly, for any N m and iV e there is a negative meson mass-squared direction. However, 
note that its value is of lower order in 1/Mp than the leading term (which is of zeroth 
order) in the negative squark mass-squared that is driving the quark phase transition; 
see (15.131) . Hence even without calculating the critical temperature T]f, that would cor- 
respond to m 2 , we can be sure that it would be much smaller than Tf? of Subection 15.21 
Therefore, by the time the temperature starts approaching Tjf, the system would have 
already undergone the second order phase transition in the quark direction and would be 
rolling along the Q axis. So let us get back to considering (16. 5p with Q ^ 0. 
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One easily finds that, up to order 1/Mp, the two eigenvalues of the matrix M are: 

1 



m 2 , = 2h 2 Q 2 N 2 m + O 



M 2 



ml = k2(Q2 ^7 NmNe + ^1^2 W + ^ 2 Q\Q 2 ~ » 2 ?Nl 

+ Ah 2 fi 4 Q 4 N m N e -h 2 fi A (Q 2 + fi 2 ) 2 N e ] . (6.9) 

The expression for m 2 _ seems divergent for Q = 0. However, we just saw in the previous 
paragraph that at the point Q = things are completely regular.^ The reason for 
the apparent problem in ( 16. 9j) is the following. To obtain the last formulae, we had to 
expand the square root in (16.51) for small 1/Mp. This is perfectly fine when Q is of order 
[L. However, when Q — > one has to be more careful as the "leading" zeroth order 
contribution, 2h 2 Q 2 , in the above mass formulae becomes <<c /j, /Mp (keep in mind that 
/j, 2 /Mp may be small, but it is definitely finite). So to extract the correct behaviour of 
( 16.51) in the limit Q — > 0, one has to first expand in small Q and then expand this result 
in small 1/Mp to the desired order. Doing that, one recovers (16.81) at the zeroth order in 
the Q-expansion. 

The lesson we learn from the above considerations is that (16.91) is valid in a neighbor- 
hood of the point Q = fi, in which Q and fi have comparable orders of magnitude. In 
that neighborhood both m\ are positive definite except at Q = fi, where one hasH 

9 ,, 9 2W 2 N e 2h 2 u 6 N 2 (N m -l) 

So, in principle, the sign of m 2 _|g =At depends on the relative magnitudes of W and hfi 3 , 
except for the case N m = 1. However, if one wants the positive vacuum energy density 
in the metastable state to be very small, then the following relation should be satisfied 

mini: 

W 2 

||^yiV e , (6.11) 

where ~ means that the two sides are of the same order of magnitude. If one assumes this 
relation, then the first term in m 2 _ is dominant as by order of magnitude it is ~ h 2 fi 4 /Mp 
and so m 2 _\Q =il < 0. Hence, when (16. lip holds, the new effect due to Wo ^ corrects the 
order 1/Mp results instead of those at 0(1/Mp). 



30 This is, in fact, the main reason we performed that computation explicitly; it was already clear from 

()6.ip that if there is a critical temperature at Q — 0, then it must be that Tjf <<c . 

31 Strictly speaking, (|6. 10[) is also valid in the small neighborhood in which Q — \i — > 0, i.e. when 
Q — [i <K n 2 /Mp. The argument is analogous to the one below eq. (|6.9p . 
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Let us see what conclusions one can make for the sign of mt throughout the interval 
[0, n] for the cases when the Wq contribution is of order 1/Mp and of order 1/Mj>, respec- 
tively. In_either case, the sign of the eigenvalue m 2 _ is determined by the sign of detM, 
see (16.5 



2j However, in the first case, the leading order contribution is ~ Q 2 {Q 2 — A* 2 ) 2 
and so is positive-definite except at Q = and at Q — fj, At these two points the sub- 
leading term (of order 1/Mp) determines the sign and it is negative due to the negative 
Wq contribution. On the other hand, in the case when ( 16. lip holds we should only look 
at the terms of up to 0(1/Mp). The leading contribution to detM, which is first order in 
1/Mp, is now detM« ~ Q 2 (Q 2 - /j 2 ) 2 h 2 N m N e - Q 2 W^N e /Mp, which upon using (I6TTT1) 
gives 

detM« « Q 2 h 2 [(Q 2 - n 2 ) 2 N m N e - /i 4 iV 2 ] < (6.12) 

at any point in (0,/i] as N m < N e . To recapitulate: when (16. lip is imposed, we find that 
m 2 _ < for any Q G [0, /x], while without (16.111) vn?_ is negative only in small neighborhoods 
of Q = and Q = fi. In a similar way one can show that m+ is positive-definite without 
(16. lip , whereas with (16. lip it is negative only in a small neighborhood of the origin. 

The above conclusions imply that at some tempreature Tf , below of Section 5, 
there will be another second order phase transition, this time in the meson direction 
corresponding to m 2 _. We will explain in Section [8] that this phase transition is actually 
necessary in order for the system to roll towards the metastable vacuum, due to the shifting 
of the latter from (ip 1 ,ip 2 ) = (0,0) to ((p 1 ,<p 2 ) ~ (^/Mp, fi 2 /M P ). However, computing 
the critical temperature is now much more complicated than in the quark case. One 
reason is that the phase transition occurs at some Q ^ 0. (With (I6.1ip imposed, this could 
apriori be any point in the interval [0,/i], whereas without this constraint it has to occur 
at Q = /!.) Another, much more serious issue is that for temperatures below there are 
masses that are greater than the temperature and so the high-temperature approximation 
(13.21) cannot be used. Instead, one should consider the full integral expression: 

Veffix) = V (x) + ^J2 ±n i dx x " ln ( X T exp(- v /x 2 + mKx)/T 2 )^ , (6.13) 

where nj are the numbers of degrees of freedom and the upper (lower) sign is for bosons 
(fermions). 

Unfortunately, that means that we cannot obtain a simple analytic answer for the 
meson critical temperature T£. However, we can make an estimate of its magnitude in 



32 If detM < 0, then the expression under the square-root is greater than (TrM ) 2 and so m 2 _ < 0. 
33 More precisely, except in very small neighborhoods of those points; see the remark in footnote I3T1 
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the vein of [TJ. Namely, let us consider a path in field space connecting the point along 
the Q-axis, at which the rolling in the meson direction starts, with the point (Q, (fx, ^2) 7^ 
(0, 0, 0), where the metastable minimum is, and take into account only fields whose masses 
change significantly along this patho Now, let us assume that all relevant masses are 
either much smaller or much greater than T. We will see below that our results are 
consistent with this assumption. Then for fields with m « T we can still use ( 13. 2p . 
whereas for fields with m » T we can utilize the approximation (see (3.5) of [TJ): 

r°° / \ / m \ 3/2 / m \ 

± dxx 2 \n{l =f exp(- v / x 2 + m 2 /T 2 )J ~ T 4 [YtJ exp \t) ' ^ 6 ' 14 ^ 

In principle, the full temperature-dependent part of the effective potential is obtained by 
summing over all fields. Note however, that the exponential in (I6.14p strongly suppresses 
the contributions of the fields with mass m ys> T. In other words, as long as there is at 
least one field with m «T, one can neglect the heavy fields to leading order. 

As we saw above, without imposing (16.111) the phase transition occurs at Q = \x. Let 
us assume that this is true also with (16.111) . Clearly, for the latter case this assumption will 
only give us a lower bound on Tjf, but this is as good as one can get in that case without 
studying the full dynamical evolution of the system. At Q = fi, the heavy masses in our 
system are raf^ , rn 2 m , m 2 Req2 _, mj m and m 2 + . All of them are ~ h 2 ^ 2 + 0{1/Mp) 
and to leading order remain constant along the above field space path. The light masses 
that determine in this approximation are m 2 Reqi _, m 2 mqi _, m|j e(?2+ , m 2 Imq2 _ and m 2 _ IlI 
Hence we can obtain an estimate for the critical temperature by solving 

d%y l T = -m 2 _ , d v _V}. = -d v _V at (Q, <p u ip 2 ) = (ji, <p{, ip* 2 ) . (6.15) 

In this system of equations denotes the temperature-dependent part of the effective 
potential, that results from (13. 2p by summing only over the light fields, and is the 
linear combination of tp\ and <y9 2 that corresponds to the mass-squared eigenvalue m 2 _. 
One then finds that 

~ w V— (6 16) 

c hM 2 M P ' 1 J 

where in the second equality we have used ( 16.111) . Note that the heavy masses (let us 
denote them by with h running over all heavy fields) are all of zeroth order in 1/Mp, 



34 Recall that the reason for this is that fields with (nearly) constant masses contribute only to the 
field-independent T 4 term in the effective potential, which we are not interested in. 

35 All of them vary significantly (meaning that the magnitude of their change is comparable to the 

magnitude of their leading order) along the path of interest and so contribute essentially to the effective 
potential. 
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whereas the light masses (denote them by mi with I running over the light fields) are 
all multiplied by the small constant h compared to (16.161) . So our assumption, that at 
temperatures near Tf we have rrih y> T >>> mi, is consistent with the estimate in ( 16.161) . 

Before concluding this section, let us comment on the first order phase transition found 
in j8]. Its critical temperature is of order (T^) 2 ~ hp 2 . When it is reached, tunnelling 
becomes possible between ($) = and a minimum away from the origin in the meson 
direction ($) ~ Ijv., which at zero temperature becomes the supersymmetric vacuum at 
(<&) = h~ 1 p e i 3N ~ N f)/( N f- N ) H N On the other hand, our second order phase transition 
occurs in the field direction (p_ = ip 2 + -^-^rfi + ^(jjt) an d a t (Qi) — (Qi) 7^ ® — 
Which one of and is greater depends on the relative magnitudes of h and (p/Mp) 2 . 
However, regardless of that, the second order phase transition is much more likely to take 
precedence since the first order one, as any tunnelling event, is exponentially suppressed. 
To gain a better understanding of the phase structure in the meson direction and, in 
particular, to be able to estimate the supergravity corrections to the height of the potential 
barrier relevant for the first order phase transition, we would need to take into account 
the non-perturbative dynamically generated contribution to the superpotential. We leave 
that for future research. 

7 Towards KKLT-ISS at finite T 

The proposal of [10] is a significant progress towards finding dS vacua in string theory 
with all moduli stabilized. However, the uplifting of their AdS minimum to a de Sitter 
one has been rather difficult to implement in a controlled way. It was shown recently in 
[Tij [T2] . that this can be achieved easily by using the ISS model as the uplifting sector. 
They considered the following coupling: 

W = W ISS + W KKLT , K = K ISS + K KKLT , (7.1) 

where Wiss an d Kiss are given by (12. ip and (12.21) . whereas: 

Wkklt = W Q + ae~ bp and K KKLT = -3 ln(p + p) . (7.2) 

In the string context, the constant Wq is due to nonzero background fluxes. By tuning it 
suitably, one can achieve an almost vanishing cosmological constant and a light gravitino 
mass [IT] , which is an important improvement compared to the models with D-term 

36 This expression for </?_ is valid with or without (|6.11[) as Wo appears for first time at order 1/M%. 
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uplifting. Therefore, it would be quite interesting to investigate the phase structure of 
the KKLT-ISS model (17.11) at finite temperature. In this section, we limit ourselves to a 
discussion of the fate of the (Q, (fx, (p%) = (0, 0, 0) ISS minimum when the KKLT field p 
is included. 

In Subsection 14.11 we noted that the presence of terms linear in the meson vevs (pi and 
ifi2 can shift the minimum of the effective potential away from the origin of the ISS field 
space. This, of course, refers to terms that are not multiplying Q 2 since Q = is always 
a local minimum as long as p <<C Mp, see eq. (14.101) . However, in that section such linear 
terms were not appearing. Now the situation is very different, as the inclusion of the 
KKLT sector introduces many terms linear in (fx and (p2- One source of them comes from 
the constant piece Wq in the KKLT superpotential, as one can easily convince oneself by 
looking at the form of (G I Gj) in (I3.2ip ] 37 l In addition, there are many more contributions, 
linear in the meson vevs, from mixed terms between the KKLT and ISS sectors in the 
total TrM 2 . For example, the total fermionic TrM? is now 



TrM 2 = (e G \K AB K^ u (y A G c + G A G C )(V B G D + G B G D ) - 2j ) , (7.3) 

where A = {p, 1} and the index / runs over the ISS fields. Writing this out, we have 

TrM) = {e G [KPPK^{W p G p + G p G p ){V- p G- p + G- p G p ) 
+ 2K^K IJ (V P G I + G p Gi) (V- p Gj + G- p Gj) 
+ {G IJ + G I Gj){G IJ + G I G J )-2}), (7.4) 

£3 The last line is the familiar 



where in the second line we have used that V p (jj = V iG pi 
ISS plus supergravity contribution, which for W = Wiss + Wkklt gives the linear term 
mentioned in the beginning of this paragraph. Let us now take a more careful look at the 
remaining terms. One can verify that the second line contains, among many others, the 
terms 

( t -^[WKjW 1 + WK'Wj] ) . (7.5) 

[p + py 

It is easy to see that these also lead to contributions which are linear in tp\ and ip2- Further 
linear terms come from the first line of (17.41) and also from the mixed terms in the scalar 
potential 

V = e K (K AB D A WD B W -3\W\ 2 ) (7.6) 



37 Recall that the full expression is (e G G I Gi). Hence the mixed terms in (|3.21[) give contributions 
proportional to ifiWo and to ip% Wq. 

/inn \- r-\ s~i - I , , 

P — U P 

,>p an< i r^ p 



38 The relation V P G/ = V/G p is due to <9/G p = d p Gi = - W p Wi/W 2 together with the fact that the 
only nonvanishing Christoffel symbols are Y p and T p 
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which determines the bosonic TrM^ via [22] : 



8PV 



TrM fe 2 = 2 (K CD — '^j-g ) . (7.7) 

Therefore, in the coupled KKLT-ISS model, generically, the origin of the ISS field 
space (Q,(fii,(f2) = (0,0,0) is not a local minimum anymore. The coefficients of the 
terms linear in the meson vevs are functions of the KKLT field p. Hence (p) is related to 
the magnitude of the shift of the minimum in the ISS plane. Although very interesting, 
we leave further analysis of the KKLT-ISS system for future work. 



8 Discussion 

We studied the effective potential at finite temperature for the ISS model coupled to su- 
pergravity. Assuming that at high temperature the fields are at the origin of field space, 
which is a local minimum, we investigated the phase structure of the system as it cools 
down. In the quark direction, the situation is analogous to the rigid case [8]. Namely, 
there is a second order phase transition at certain critical temperature, T®. The effect 
of the supergravity corrections is to decrease compared to its global supersymmetry 
counterpart. In the meson branch however, a new feature appears whenever the super- 
potential contains a nonvanishing constant piece Wq. Recall that in the global theory all 
meson fields always had a local minimum at the origin of the meson direction, with no 
tree- level contributions to their masses-squared. Now we find that, when Wq ^ 0, for some 
of them this ceases to be true at some temperature Tjf, below which negative tree-level 
supergravity corrections to their effective masses-squared are outweighing the positive 
one-loop temperature dependent contributions. Hence, the supergravity interactions lead 
to the occurrence of a new second order phase transition whenever Wq ^ 0. 

Since T£ < T®, as we saw in Sectional the fields first start rolling away from the origin 
in the quark direction. When the temperature decreases enough, the same happens also in 
the meson direction. However, unlike in the rigid case, the second phase transition does 
not imply that the system is moving away from the supersymmetry breaking vacuum. 
The reason is that the coupling to supergravity leads to slight shifting of the position 
of the metastable minimum [12] O Whereas in the rigid theory it was given by (gi) = 
(qi) = ^ I-Mm an d ($) — 0, in the locally supersymmetric case some of the meson vevs 



39 This is not a trivial consequence of including supergravity, as in the ISS model coupled to sugra with 
Wo = there is no shifting of the metastable vacuum compared to the global case. 
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also acquire nonzero value: ((diag </>n) , {diag 4>22)) = (fi,<f2) ~ (n 2 /Mp, fi 2 /Mp) <K (qi). 
Hence the latter phase transition is, in fact, a necessary condition for the system to evolve 
towards the metastable vacuum. Of course, to follow with more precision the evolution 
of this system as the temperature decreases, one has to study the full effective potential 
for both (q) 7^ and ($) 7^ away from a small neighborhood of the origin. 

In the above paragraph, we reached the conclusion that the final state of the system at 
T ~ is likely to be the metastable vacuum. However, one should be cautious since, sim- 
ilarly to [8], our considerations assume thermal equilibrium. Hence, although suggestive, 
they are not completely conclusive. As was pointed out in [8], for proper understanding 
of the evolution of the system one should address also the dynamics of the fields at fi- 
nite temperature. Actually, even before worrying about dynamics, one may be concerned 
that in our case the situation is complicated by the existence of new supergravity-induced 
super symmetric minima@ Indeed, it was shown in pj] that such a minimum occurs in 
the KKLT-ISS model@ However, this new vacuum only appears due to the interaction 
with the KKLT sector; it is easy to see that the last condition in the solution for this 
minimum, eq. (19) of [TT], is only satisfied with an appropriate choice of (some of) the 
tunable KKLT parameters Wq, a and b. Still, one may wonder whether there could be 
a solution to the supersymmetry preserving equations for the ISS plus supergravity sec- 
tor alone. In Appendix B we show that this is not possible (at least in the field-space 
drections of interest, i.e for an ansatz for the vevs that is of the same type as the one in 

CI]). 

The new solution of [11] is only one indication that the KKLT-ISS system is quite 
intricate to study. Another is the fact that, as we saw in Section [TJ the interaction 
with the KKLT sector leads to shifting of the high-temperature minimum of the effective 
potential away from the origin of the ISS field space. Understanding the phase structure 
of this system is of great interest. However, the technical complications involved are 
rather significant. Therefore, it may be beneficial to gain preliminary intuition about it 
by considering the recently proposed O'KKLT model [30], as the latter is much simpler 
while having quite similar behavior. In addition, the O'KKLT model has a significance of 



40 That is, susy minima other than those induced by the non-perturbative superpotential Wads, see 
eq. (|2.8[) . Recall that the Witten index gives only the number of global susy minima of globally super- 
symmetric theories. Hence, in the present context it is not applicable and so one cannot immediately 
rule out the presence of additional solutions. 

It is true that this minimum is much further out in field space than the metastable one, but its very 
existence raises the possibility that it may be quite premature to make conclusions about the final state 
at low temperature, based solely on studying the immediate neighborhood of the origin. 
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its own, as it was argued in [30] to be of value in studies of cosmological inflationo We 
hope to address this in the future. 
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A Useful mass matrix formulae 

By differentiating the F-term part of (14.11) w.r.t. to x 1 an d Xj one finds: 

d T d J V = e G [M P {G LJ + G L G J )(G LI + G L G I )-M P G I G J + 5j(M P G L G L -2M P )} . (A.l) 

Substituting e G = e K/M p\W\ 2 /M P and 

r Mr r 1 W LI + K L Wj + KjW L K L Kj 

we see that all powers of W in the denominator cancel out. We want to compute the value 
of the resulting expression for two cases: One is for zero background vevs of all scalars. 
And the other is for cpi, <f 2 — but Q ^ 0. 

In the first case, i.e. for Q, tpi, if 2 = 0, one has that (K) = and (Kj) = for V/. Let 
us denote (W) = Wo and for the moment consider Wo 7^ for more generality. We find 
that 

(djd J V) = (W IL W JL ) + [5j{W L W L ) - (^/W J )) , (A.3) 

where all higher orders in the 1/Mp expansion vanish due to (K), (Kj) = regardless 
of the value of Wo- Realizing that the vevs of all double derivatives of W vanish at the 
origin of field space, we finally obtain: 

(djd J V) = -L ^(WjV 1 ) - (WjW J )) . (A.4) 

Note that this result, apart from being exact to all orders, is also completely independent 
of W - 



42 In fact, this statement only applies to a variant of O'KKLT considered in Section 3 there. This same 
variant is also the model that is most useful as an approximation of KKLT-ISS since it is the one, whose 
supersymmctric global minimum is at finite field values. 
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Applying (1A.4j) to compute the mass matrices m 2 = , where q is either q\ or q\ and q 



is either q\ or qi, one arrives at half of the relations in (I5.12p . The other half, i.e. of 
types m 2 qq and m? 5 , can be derived in a similar way. (It is perhaps more convenient to 
differentiate the potential in the form (13.41) .) One finds that at the origin of field space: 

(d K d L V) = (K IJ W KLI Wj) , (A.5) 

which is exactly the expression for the rigid case, since all second derivatives of the 
superpotential vanish for zero background fields. And again, this is exact result to all 
orders in 1/Mp. 

In the case of ipi, (f2 = but Q ^ 0, we will be interested in the meson mass matrices. 
For /, J running over the meson field components only, we still have that (Kj) = 0, 
although now (K) ^ 0. Keeping W ^ and expanding (lA.lj) up to order 1/Mp, we find: 



(djd'V) = (W IL W JL ) (A.6) 



+ — f Uj(W L W L ) - (WW 3 ) + (WjK L W JL + W IL K L W J + KW IL W JL )\ 

IVl p v ) 

l ¥ ^kwJv 1 ) -2W^+ (k)(WjK l W jl +w il k l W j ) + {^-w IL w JI ^. 



+ 



Obviously, for Q = the above formula agrees with flA.4j) . However, recall that the /, J 



indices in it run only over the meson fields, whereas L runs also over the squarks. 

Finally, one can easily convince oneself that if K, L run only over the mesons, whereas 
the squarks are the only fields with nonzero vevs, then (c^-c^V^) = to all orders and for 
any value of Wo- 



B On existence of new susy solutions 

We investigate here whether there are solutions to the supersymmetry preserving equa- 
tions for the ISS model coupled to supergravity, with the non-perturbative superpotential 
Wads of eq. (I2.8P still neglected. 
The susy equations are: 

D qi W = 0, D qi W = 0, D^ u W = 0, D 022 W=O. (B.l) 

Similarly to [llj, we consider the following ansatz for the expectation values of the quark 
and meson fields in the solutions we are seeking: 

(qi) = fiil Nm , (qi) = /i 2 l7V m , ($n> = Vi%N m , ($22) = v 2 \ Ne , (B.2) 



33 



where all vevs are real. Evaluating (IB.lj) in this background gives: 



D qi : hfjki/! + in(W) =0 
D 011 : /x(/xi/x 2 - xx 2 ) + Vi(W) = 



% : h/nui + n 2 (W) = 



-V + ^W = 0. (B.3) 



The two equations on the first line of flB.3|) imply that /x 2 = /x 2 . On the other hand, the 
ones on the second line lead to: 

,2 



/il/i 2 - /! 

Now, using this last relation and the equations for D qi and -D$ 22 , we find 

1/2 



(B.4) 



/ii 

M2 



(/Xi/X 2 ~ fl*) , I/ 2 



/'I 



(B.5) 



So far, we have expressed all vevs in terms of one of them, which could be either /xi or 
/x 2 . Let us choose this to be To obtain an independent equation for it, we need to use 
the explicit vev of the superpotential: 



(w) = h (niiwiN m - n 2 (yiN m + ^N e )) . 

Combining this with ( IB. 41) and the D qi equation in ( IB. 31) . we find: 

H 2 N e 



(B.6) 



M2 



Mp 



ViH2N m — /x [ N, 



(/Xi/x 2 - fi 2 ) 



(B.7) 



where we have reinserted the explicit dependence on Mp that comes from DjW = diW + 
(Kj/Mp)W. Let us now consider first the case /x 2 = /Xi. Then (IB. 7ft becomes a quadratic 
equation for /x 2 , whose solutions are: 



Mr 



Mi 



2N n 



2N n 



1-4 



Mi 



N N 



(B.8) 



Since we would like both /xi <3C Mp and /x « Mp in order to have a reliable field theory 
description, only the plus sign in (IB.8I) is meaningful. Hence, we have 



Mr 



// 



Mi 



// 



Mi 



N e + 



Ml 



(B.9) 



43 Note that, up to now, we have not used the explicit form of (W). So our results (|B.5[) . together with 
jtxf = xxf , are valid also for the KKLT-ISS set-up considered in [11] . In fact, the solution given in their 
cq. (19) is only valid for /xi = /X2 , in which case one can see that it agrees with (|B.5[) . 
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This implies that fif — fi 2 < 0, which is inconsistent with (IB. 51) since the vevs 1/1 and 
are real. If we take in turn \ii = —fii and repeat the above steps, we end up with 



> 2 / 4 

/Xl \ / /i \ /nil 1 



+ Kti , (B.10) 



which is again inconsistent for real vevs. So we conclude that coupling to supergravity 
does not increase the number of vacua of the ISS model. 
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